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1. INTRODUCTION 

The increasing growth of load demand, lack of active power generation as well as technical 
constraints in power system are the challenges faced by power utilities in providing reliable power to 
consumers. Increase in load demands from the industrial, commercial and residential sectors has caused 
congestion in power network and this situation leads to voltage instability in power system [1]. Unstable 
power system due to voltage decline may lead to voltage collapse. Therefore, maintaining voltage stability is 
a major concern in power system operation. An approach to improve voltage stability in power system is by 
integrating renewable energy based distributed generations (DGs) which are usually installed in power 
distribution system [2-5]. However, the integration of renewable energy based DGs into a power system 
causes certain technical effects because the network design initially does not consider the integration of DG. 
The rapid growth of DG installations and high penetration of DG can cause many technical issues such as 
power quality. Integrating DG at random locations may cause increased power losses and voltage instability]. 
To maintain voltage stability and reduce power loss, it is important to place DGs optimally in a power 
system. 

Several methods have been developed for determining optimal DG placement using either analytic 
or heuristic optimization approaches. In [6], power stability index is used to determine the critical buses in a 
system where these buses are considered as optimal DG locations. A modified voltage index is used as an 
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objective function to improve voltage stability margin and mixed-integer nonlinear programming is applied 
to solve the optimization problem [7]. In addition, the voltage stability index (VSI) has been developed to 
determine the optimal placement and sizing of DG in a radial distribution system considering load growth 
[8]. Among the applied heuristic optimization techniques, the particle swarm optimization (PSO) has been 
used for determining the size and location of multi DGs considering multi-objective index [9]. In addition, 
considering heuristic optimization, the cuckoo search algorithm [10] and a constructive heuristic algorithm 
[11] have been used for solving the optimal DG placement problem. In all the above-mentioned methods, a 
general DG model with real power injection has been considered without considering any renewable energy. 

One type of DG that is widely used in modern power system is the photovoltaic based distributed 
generation (PVDG). The PVDG system provides additional active power into the network, thereby 
facilitating reverse power flow and reducing energy consumption. The reactive power is solely provided by 
the conventional generators and therefore, reactive power limitation is neglected when finding the optimal 
PVDG location and sizing. To improve voltage stability margin and reducing power loss in a power system 
with PVDGs, appropriate planning of PVDGs is considered by optimal placement and sizing of PV. Some 
research works have been done for determining optimal placement and sizing of PVDG considering voltage 
stability improvement [12-14]. In [12], the PSO technique was used to determine optimum placement of 
PVDG and wind turbine by incorporating voltage stability factor. In [13], the artificial bee-colony 
optimization technique is developed for determining optimal sizing of hybrid PV/wind turbine/fuel cell 
system by maximizing the voltage collapse point using the VSI. In [14], the self-correction algorithm was 
used in determining optimal sizing of multiple PV and battery storage units considering voltage stability 
margin. 

In this paper, a relatively new heuristic optimization technique known as wind driven optimization 
(WDO) [15] is applied for solving the optimal placement of PVDG units. In the proposed optimization 
problem formulation, the objective function is to maximize the voltage stability margin. Voltage stability of a 
power system is assessed by using a relatively new voltage stability index known as the improved voltage 
stability index ([VSI) [16]. Thus, for determining the optimal location of PVDG, the objective function 
considers maximizing the IVSI. 


2. THE MATERIAL AND METHOD 

This section describes [VSI for voltage stability analysis, PVDG modelling, and application of 
WDO for solving the optimization problem. The IVSI is an indicator for assessing whether a power system is 
stable from the voltage stability viewpoint. It is used in the objective function in finding the optimal location 
of PVDG for improving voltage stability in a power system. This section also discusses on the application of 
WDO technique in determining PVDG optimal placement. In this study, the entire modeling and simulation 
is carried out using MATLAB software. 


2.1. IVSI for voltage stability assessment 

The development of [VSI for voltage stability assessment is based on the current and power flow at 
a bus in a power network. The mathematical formula for [VSI for n-bus system is derived in [16] and it is 
given by, 
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where P; + Q; is the load active and reactive powers at the receiving end bus, V;is the voltage at the sending 
end bus, 6, is the phase angle difference between sending end and receiving end buses and Gj and Bj are the 
elements in the bus admittance matrix. 

The IVSI is an indicator to assess voltage stability at the buses in a power system. The IVSI has a 
value between 0 and 1. A bus is said to be stable if the IVSI value is less than 1 (IVSI <1). However, a bus is 
said to experience a voltage collapse if the [VSI value equals to 1. When running a simulation to assess 
voltage instability using IVSI, the loading margin at a bus is increased gradually until the [VSI value reaches 
approximately 1. 
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2.2. IVSI for voltage stability assessment 
In this study, the static model of PVDG is based on a single-diode equivalent circuit [12] where the 
active power is expressed as, 
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where G is the solar irradiance set at constant 1000 W/m?, Tc is the solar operating temperature, Voc is the 
open circuit voltage, sc is the short circuit current, Npy is PV module number, TA is the ambient temperature 
which sat at constant 35°C, NOCT is the nominal operating cell temperature, K; is the short circuit current 
coefficient, K, is the open circuit voltage coefficient, 7,7 is the reference operating temperature at 25°C. All 
the electrical characteristics at standard test condition (STC) are taken from the PV datasheet provided by the 
manufacturer. 


2.3. Formulation of the Optimization Problem 

The optimization problem formulation considers maximizing voltage stability margin of a power 
system. The search for optimal solution is subjected to the following objective function and technical 
constraints. To determine the optimal PVDG location, the weakest bus in a test system has to be determined 
by calculating [VSI at each load bus. From the mathematical formula in (7), the [VSI value is between 0 and 
1, such that if a bus has an IVSI value of 1, it means that the bus is the weakest bus from voltage instability 
viewpoint. Thus, the [VSI at each load bus as an objective function need to be maximized for determining the 
weakest bus. The mathematical formula of objective function for optimal PVDG location is given by, 


f, =max {IVSI, ' (6) 


The technical constraints considered in the optimization problem are the PVDG, voltage and power 
flow constraints. The PVDG is limited by the available power sources at any given location. Thus, it is 
necessary to constraint the capacity between its maximum and minimum level of active power, which is 
given by, 


| ES od (7) 
The system operation should be kept within +5% of voltage deviation and the voltage constraint is given by, 

0.95 pu <V, <1.05 pu (8) 
The power flow equations to be satisfied are given by, 
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where JN is the number of bus, Pipvng and Q;pvng are the PV generated real and reactive powers at bus i, 
respectively and P;,p and Qj pare the real and reactive power demands at bus /, respectively. 


2.4. Application of WDO for Optimal Location of PVDG 

WDO is based on the Newton’s second low of motion which is used to describe the motion of air 
parcels on the earth’s atmosphere [15]. In the WDO topology, the position of small air parcels with random 
velocities in an N-dimensional space is updated based on the physical equations that govern large-scale 
atmospheric motion. The movement of an air parcel in one dimension is contributed from the movement of 
air parcels in another dimension. The velocity of air parcels is derived from the concept of an air parcel 
moving with the wind in which the updating velocity is given by, 
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where Uney and Ugia are the updated velocity and current velocity, respectively, Xma, is the distance from the 
highest-pressure point, Pinay is the maximum pressure, P.ig is the pressure at the current location, T is the 
temperature, g is the gravitational acceleration, and R, a, and c are constants. 

The third term of equation (13) is a force towards the location of maximum pressure which is 
represented by the global best location of the optimization problem. Once the new velocity is calculated, the 
position of air parcel will be updated using the following equation. 


Xnew = Xoa + (new X At) (14) 
where Xoi¢ and Xnew are the current position and new position of the air parcel, respectively and the time step 
At is assumed as 1. 

The WDO is applied for solving optimal location of PVDGs and it is described in terms of a 
flowchart as shown in Figure 1. From the flowchart, the algorithm starts with the initialization stage where all 
input parameters are defined to establish the problem boundaries. After setting the input, the first step of the 
optimization problem is to find the optimal PVDG location. In this algorithm, the PVDG is represented as an 
air parcel. During the simulation, the pressure of each air parcel (objective function) is evaluated at its current 
position. Then, the population is ranked based on its pressure and the velocity is evaluated using equation 
(13). The position of the air parcel for the next iteration is updated using (14). This procedure continues until 
it reaches the maximum number of iteration. Finally, the best location of the last iteration is recorded as the 
optimization result. 


2.5. Application of PSO for Optimal Location of PYDG 


The PSO is an evolutionary computation algorithm based on a social-psychological metaphor [17]. 
In PSO topology, the best individual position of particle, i is represented by, 


Pbest, =(Pbest,, + Pbest,, + Pbest,, +...+ Pbest,, ) (5) 
whereas the best position of the entire swarm in n-dimensional search space is represented by 
Gbest, = (Gbest,, + Gbest,, + Gbest,, + ...+ Gbest,, ) (16) 


An individual in a swarm approaches the optimum through it present position, present velocity, 
previous experiences of Pbest and Gbest and the experience of its neighbours. The particle of new position 
and velocity are determined using the following equations, 
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Where i = 1,2,3,...,1 is the number of dimension of each particle, c; and cz are weight factors, r; and rz are 
random number between (0,1) and k is the number of iteration times. The inertia weighting factor is 
represented by w of the following equation, 


omax—o@min . 
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Figure 1. Implementation of WDO algorithm for 
solving optimal location of PVDG 
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Figure 2. Implementation of PSO algorithm for 
solving optimal location of PVDG 


The PSO technique is applied for solving optimal location of PVDGs and it is described in terms of 
a flowchart as shown in Figure 2. The first step of the optimization problem is to find the optimal PVDG 
location. In the PSO algorithm, the PVDG is represented as a particle in the swarm population. The algorithm 
starts by initializing all the input parameters as shown in the flowchart. During the simulation, Popes: (objective 
function) is evaluated at its current position using equation (15) whereas Goes; is evaluated using equation 
(16). The new position and velocity of particle are updated using equation (17) and (18). This procedure 
continues until it reaches the maximum number of iteration. Finally, Goes; of the last iteration is recorded as 
the optimization result. 
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3. RESULTS AND DISCUSSION 

The simulation results for optimal placement of PVDGs in the IEEE 30 bus, IEEE 118 bus and 69- 
radial distribution systems using WDO and PSO are presented and discussed. Simulations were carried out to 
determine the weak buses in the test systems using IVSI. During the simulation, each load bus is evaluated 
individually by maximizing the IVSI as the objective function. The load factor is increased by a certain 
percentage of the base value until IVSI approaches 1, which is the voltage collapse point. The first bus that 
experiences voltage collapse is said to be the weakest bus in the system and it is marked as an optimal PVDG 
location. Thus, the PVDG is installed at this bus. For multiple numbers of PVDG, the first PVDG will be 
integrated using the bus ranking approach. The first PVDG will be integrated at the weakest bus (e.g. Bus 61 
for the 69-radial distribution system) whereas the second PVDG will be integrated at the second weakest bus 
(e.g. Bus 64 for the 69-radial distribution system). Similar procedure is applied for the IEEE 30 bus system 
and IEEE 118 bus system. 

Table 1, Table 2 and Table 3 show the IVSI values with variable load factors for the 69-radial 
distribution system, IEEE 30 bus, and IEEE 118 bus, respectively. The load factor of each bus is increased in 
percentage from the base load value. Results in Table 1 show that the optimal PVDG location is Bus 61 
because at load factor of 190%, the IVSI value equals 1.0 and Bus 61 is said to experience a voltage collapse. 
However, other buses in the system still operate within voltage stability limit with IVSI values less than 1. 

Results in Table 2 show that the optimal PVDG location of the IEEE 30 bus system is Bus 30. It is 
the first bus to reach voltage collapse point (IVSI=1) at a load factor of 150%. The next optimal PVDG 
location is Bus 26 and then followed by Bus 5. 


Table 1. IVSI with load increment for the 69-radial distribution system 








WDO PSO 

ET EAGER) Bus 65 Bus 64 _—_— Bus 61 Bus 65 Bus 64 Bus 61 
100 0.2622 0.4092 0.4316 0.2622 0.4092 0.4316 
120 0.3286 0.5127 0.5400 0.3286 0.5127 0.5400 
140 0.4079 0.633 0.6635 0.4018 0.6266 0.6591 
160 0.5211 0.7714 0.8003 0.4833 0.7534 ~=0.7910 
180 0.5936 0.9100 0.9718 0.5750 0.8960 0.9389 
190 0.7093 0.9842 1.0 0.6256 0.9746 1.0 





Table 2. [VSI with load increment for the IEEE 30 bus system 








WDO PSO 

TAL EBCEOE C7) Bus5  Bus26  Bus30  Bus5  Bus26 Bus 30 
100 0.3682 0.1600 0.2057 0.3500. 0.1479 0.1959 
120 0.1801 0.3057 0.3191 0.4417 0.2371 ~—-0.3140 
140 0.5984 0.4473 0.5894 0.5589 0.4237. «0.5772 
150 0.6703 0.7019 1.0 0.6583 0.6958 1.0 





Results of the IEEE 118 bus system depicted Table 3 show that the first bus to experience voltage 
collapse is Bus 117 and then followed by Bus 1 and Bus 13. Thus, Bus 117 is considered as the optimal 
location for PVDG. From the results, each test system gives similar optimal PVDG location using the two 
different optimization techniques; WDO and PSO. However, the [VSI value at each load factor varies with 
different optimization techniques. Evaluation for the best optimization technique will not be performed 
considering that the differences in the IVSI values are not significant. However, the two different 
optimization techniques; WDO and PSO are compared in terms of convergence characteristic. Figures 3, 4 
and 5 show the convergence characteristic for optimum PVDG location for the 69 radial-distribution system, 
30 bus system and 118 bus system, respectively. The results in the figures show that the WDO gives faster 
convergence with minimum objective function values compared to the PSO technique. 


Table 3. IVSI with load increment for the IEEE 118 bus system 








WDO PSO 
onda) eid: = Bagh Raed JBueis Bust <“Buvat7 
100 0.0324 0.0102 0.0895 0.0324 0.0102 0.0895 
120 0.1106 0.0905 0.1407. 0.1010 0.0837. 0.1347 
140 0.1265 0.1177 0.1632 0.1196 0.1014 0.1531 
160 0.3143 0.2904 ~—-0.3879~—Ss«0.3021.—Ss«0.282—Ss«0.3733 
165 0.8379 0.9115 1.0 0.8261 0.8959 1.0 
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Figure 5. Convergence characteristic for optimum PVDG location using WDO and PSO for the IEEE 118 
bus system 


Table 4 shows the MATLAB time simulation for the implementation of WDO and PSO techniques 
in determining optimal location of FVTG generator. The WDO technique has faster time simulation 
compared to the PSO technique where differences can be seen for the IEEE 30 bus and 69 bus system 
simulations. Therefore, it can be concluded that the WDO technique is most suitable to be used in PVDG 
optimal location determination. 


Table 4. MATLAB simulation time for optimal location of PVDG 
MATLAB simulation time (s) 

IEEE 118 bus IEEE 30 bus 69 bus system 

WDO 816.703 805.418 804.531 

PSO 816.831 810.924 807.999 





Optimization Technique 








4. CONCLUSIONS 

This paper has presented the application of WDO for determining optimal PVDG location by 
maximizing voltage stability margin, IVSI as an objective function. The WDO is applied to solve the 
optimization problem and it is compared with the PSO technique. The WDO and PSO techniques are tested 
on the IEEE 30 bus system, IEEE 118 bus system and 69-radial distribution system. Comparison on the 
performance of the WDO and PSO techniques showed that the WDO technique gives faster convergence 
compared to the PSO technique. Analysis of results proved that integration of PVDG on the weakest bus 
could improve power system voltage stability by increasing the value of IVSI. Increased PVDG penetration 
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level also increases voltage stability and improves bus voltage profiles. The maximum level of PVDG 
penetration is determined by ensuring that the bus voltage does not exceed the upper voltage limit of the 
nominal voltage. 
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